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Abstract 

The Weyl semimetal (WSM) is a novel topological gapless state with promises exotic transport 
due to chiral anomaly. Recently, a family of nonmagnetic WSM candidates including TaAs, NbAs, 
NbP etc is confirmed by first principle calculation and experiments. The TaAs family are reported 
to display the large unsaturated magnetoresistance (XMR), which have not yet been explained. 
Here we give a theoretical calculation of XMR based on the extended effective-medium approach 
to Weyl semimetals. We predict the power law of XMR at high magnetic field and the ”turn-on” 
magnetic field, which are well in agreement with experiments data. Furthermore, we investigate the 
0-dependence magnetoresistance and find the transition between the postive XMR and the negative 
magnetoresistance induced by chiral anomaly, which should be confirmed by further experiments. 
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Gapless topological phases, especially the Weyl semimetal(WSM), have attracted much 
intereslfn^. Recently, TaAs-family materials have been proposed as possible candidates for 
realizing the WSM®. TaAs is a kind of nonmagnetic material with preserved time reversal 
symmetry (TRS) and broken inversion symmetry (IS), and unlike previously proposed py- 
rochlore material^, its lattice structure satishes a C4 rotation symmetry. According to the 
ab initio calculatiorP^l^ its Fermi surface is made up of 12 pairs of Weyl nodes, of which four 
pairs with tiny hole pockets (about 2 mev) are located near a point in the = 0 plane and 
the other eight pairs with electron pockets are located symmetrically in the kz ^ Q plane. 
The exotic surface state of TaAs, the Fermi arc, has also been observed experimentally by 
Arpe^MSI. 

At the same time, its reported by several experiments on the transport property 
of TaAs, NbAs, NbF^^^^that these candidate materials for WSM all show large unsaturated 
behavior of magneto-resistance (MR). To our knowledge, this near-linear behavior of MR in 
high magnetic held has not been theoretically explained up to now. 

In Fig Id, we show the experimental data of the MR of TaAs family materials from 
different group^^®^ and for comparison, we also include the MR data of WTe^^. Its easy 
to see that The MR of TaAs family is 1-2 orders of magnitude larger than that of WTe2. 
The unsaturated MR of single crystals in this family shows similar asymptotic behavior, 
namely, quadratic under low magnetic held (< IT) and near-linear under high held, with 
the transition point between IT lOT. So far, there has been two kinds of explanations for 
the TaAs familys unsaturated MR. One of them is based on the semi-classic two-band model 
with carrier compensation at n=p, and the charge resonance can explain the quadratic-type 
unsaturated behavior of MR quite well. Unfortunately, as the measurement of Hall conduc¬ 
tivity at low temperatures shows electron-type doping, i.e. n > p, MR will get saturated 
in high held and the resonance condition is not satished here, so the carrier compensation 
theory is incapable of explaining the near-linear type MR. The other explanation stems from 
Abrikosovs quantum MR modeP^, which is based on linear dispersion with all charge carri¬ 
ers occupying the lowest Landau levels. This theory is characterised by the LMRs stability 
against temperature, which, of course, cannot explain the temperature-dependent transition 
held H*. 

However, we hnd that the ehective media theory (EMT) which was hrst developed in 
disorder systems could provide a simple physical explanation of the unsaturated MR of 
WSA/pSHll. The current phase diagram for the asymptotic behavior of MR is shown in 
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Fig.le, including carrier compensation black dotted line)^^^, the Abrikosov LMR(red 
solid line)ff^, the duality theorem for 2D systems(Ohve dashed)ff^ and the classical saturated 
MR(blue short dashed)!^. Only the EMT can explain the asymptotic behavior of MR in the 
shaded area from saturated to quadratic type unsaturated. In addition, it can also give a 
quite reasonable estimate of the transition held H*. For WSM, when the magnetic held and 
electric hied are parallel, there also exists a negative MR induced by chiral anomaly, which 
has been observed in Bii_xSbj. materiald^. In this paper, we will generalize the EMT to 
systems with chiral anomaly and theoretically investigate the MRs dependence on the angle 
between the magnetic held and electric held. 

For Weyl semimetal, the total current includes two parts: the normal current jn, and the 

anomalous current ja- Thus total current j could be described by 

A f / 

j = jn + ja = (tE + —b X E+ (1) 

Where the hrst term is normal current, sigma is the conductivity, the second term is the 
anomalous Hall ehect, where b is the momentum seperation between two Weyl modes,and 
the last term is known as the chiral magnetic ehect and A/i is chiral energy shift. The 
last two terms are general anomalous current induced by chiral anomalj^. The last term 
chiral magnetic ehect is subtler and vanishes at equilibriunP^, but we could charging it 
with the external parallel magnetic held and electric helcP, and the chiral current j^h could 
be obtained by 

jch = liE ■ B)B (2) 

Where 7 = Ti is the intervalley-scattering time between two Weyl nodes, 5 f(eir)is 

the density of state at Fermi surface. The chiral energy diherence is proportional to chiral 
anomaly E-Bhy charging, and the corresponding calculation are summarised in supplemen¬ 
tary. And the total current then could be reduced in form j = a^^'^E , and the components 
of total conductivity are dehned by 

^Q!/3 ^afSrjb-ri T B^B (3) 

where b^ is the redehned anomalous Hall conductivity (which is small than a quantum 
conductivity). In general, the total conductivity does not satisfy Onsager reciprocal 
relations due to the second term violates the time-reversal symmetry. The hrst term nor¬ 
mal conductivity could be obtained by the ehective-medium approach,when the intra-node 
scattering r is considering far greater than inter-node scattering r*. 
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For Weyl semimetal candidates TaAs, eight pairs electron pockets and four pairs tiny 
hole pocket^ are simply considered as two bands model with the carrier concentration 
n = — Ylii=iPi ■ Therefore, the electron (hole) band conductivity in magnetic 

held is given by 
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Where fin, Pp are electron-pocket and hole-pocket mobility at zero magnetic held. The 
ehective medium theorj^ yields a self-consistent equation for ehective normal conductivity, 

6a,{I-T6ai) = 0 (5) 

i=n,p 

Here, 6ai = ai — a and T is an ellipsoidal depolarization tensor which depend on the 
resulting conductivity. Thus, the total conductivity is combined the normal conductivity 
and the anomalous part in Eqs.3, then the magnetoresistance MR = the 

Hall coefficient Rh = where p{H) = . 

First, we consider that H is perpendicular to E and then the total current have no the 
contribution of chiral anomaly. For WSM, the anomalous Hall conductivity is negligible due 
to \b\ << axy, then the total current equals to the normal current approximately within the 
EMT. In Figure 2a, we assume that cro,n = (^o,p, P = Pn = Pp, then n=p. The magnetoresis¬ 
tance MR is as a function of magnetic held H at difference mobilities from 5 x 
to in Figure 2a. We hnd that the power law of MR at hxed conpensation 

n=p tends to square dependence at high magnetic held, where the scaling a is 1.8 at high 
mobility and 1.3 at low mobility 5 x Furthermore, at low 

magnetic held, the scaling a = 2 when the mobility p > . Figure 2b shows 

that Rh is alway zero at hxed n=p. The asymptotic behavior of MR conhrm us the carrier 
compensation within semiclassic two band modeP^. To further investigate the carrier com¬ 
pensation within EMT, we alter the resonance condition, and show in Figure 2c the MR as 
a function of hole mobility at diherent magnetic held lOT, 30T and 60T respectively and 
settle fin = (Jo,n = o'o.p- When fip reaches ,corresponding to 

perfect compensation p/n = 1 , the MR emerges a peak at high H due to charge resonance. 
When fip is away from pn, then the MR is oh-resonance and drops down quickly. Figure 2 
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shows that the square MR of carrier compensation within two-band model is the limit value 
of EMT at high mobility and hxed n=p. 

Next, we apply the EMT to the magnetoresistance of WSM candidates. From the exper¬ 
imental estimation n > p, /i„ > /ip for TaAs familj^^, the resonance condition is dehnitely 
unsatished. For simplicity, we normalized the electron-pocket conductivity ao,n = en/i„ = 1 
and settle fin = then we study two cases: (i) at hxed fin = fip , and varying 

p/n ; (ii) at hxed p=n, and varying fin/Pp ■ We plot the calculating results with EMT 
in Figure 3. Figure 3a and 3b show the reduced resistivities pxx and pxy as a function of 
magnetic held at diherent ratios of carrier density p/n . As p/n decreasing, the scaling a at 
high magnetic held then decreases, correspondingly. Moreover, the power law of pxx is still 
square (a' ~ 2) at low magnetic held. The scaling changes from a parabolic to a less-linear 
dependence with a crossover manner characterized by a turn-on magnetic held H*. The H* 
and the scaling a at high magnetic held as a function of p/n are shown in Figure 3c. The 
turn-on magnetic varies roughly in the range from IT to lOT, which are qualitatively consis¬ 
tent with the recant experimental dataP^on TaAs, which the H* varies from 0.5T to 6T at 
diherent temperature. The corresponding power law at high magnetic held change from .5 
to 1.5 without saturation. The quasi-linear power law is around at the ratio p/n = 0.9. The 
Hall resistivity pxy in Figure 3b shows electron-type doping consistent with our assumption 
n > p. 

For case (ii) with same carriers concentration n=p. Figure 3d and 3e show pxx and pxy as a 
function of magnetic held at diherent ratio of mobility pn = Pp ■ The resistivity pxx shows 
the power law a decrease with the mobility ratio, which are the similar asymptotic behavior 
in Figure 3a. The turn-on H* and scaling as a function of the mobility ratio pn = Pp are 
shown in Figure 3f, which conhrm the same power law with variation of p/n. The slight 
diherence of in Figure 3a and 3d is that the power law at low H decreases with the descend¬ 
ing of mobility ratio, keeps no square any more. The Hall resistivity pxy in Figure 3e show 
hole-type doping characterization, even at hxed n=p. 

Finally, we investigate the angle dependence of magnetoresistance of Weyl semimetals. 
When the magnetic held H is parallel to electric held E, the the total conductivity could 
reduced to where cxo is the conductivity at H=0. Therefore, the negative 
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magnetoresistance induced by chiral anomaly is given 

MR = -( 6 ) 

CTo + 

In general case for the angle 6 between H and E, we consider the contribution of vertical 
magnetic held component H sin0 to the normal conductivity part within the EMT, and the 
contribution of the parallel component H cos 6 to the anomalous conduct ivit}^. In Figure 
4a, we show the chiral energy difference A/i between two Weyl nodes could be charged by 
chiral anomaly. Figure 4b show the magnetoresistance with normalized an,p = = 

and settle the parameter 7 = 10“^ at different angle 6 ^ = 0 to 7 r/ 2 . The 
points show the angle dependence of MR with chiral anomaly, while the dashed lines are 
referenced MR without the contribution of chiral anomaly. At small angle the MR shows 
a transition from positive MR to negative MR at low magnetic held. For the conditions 
of carrier compensation, we could give an explicit expression for the total conductivity 
~ i+p 2 sin^ e then the 0 —dependence magnetoresistance 
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MR = 


( 7 ) 


2U2 ■ 2n + COS^ 9 

Where /i' = is the effective mobility. The magnetoresistance transition could be 

obtained by = 0 . Figure 4c shows the contour plot of MR transition between PMR 
and NMR, and the yellow solid line shows the transition position as a function of angle 6 
and magnetic held H with dimensionless parameter 6 = /i'/h = 0-3- Figure 4d show the 
transition curves at diherent parameter <5 from 0.1 to 1.0. The contribution of chiral anomaly 
increases with 6 increasing, and the magnetoresistance show more NMR characterization. 
We need solve the EMT self-consistent equations for the general cases when the carrier 
compensation condition is unsatished (i.e. n > p, fin > Pp)- For experimental realization on 
TaAs family, however, the anomaly-induced negative magnetoresistance has to be carefully 
measured because it is much smaller than the giant unsaturated PMR as a angle-dependence 
background. The parameter 5 1 shows that the negative MR is not easy to detected in 

experiments. 

In conclusion, we extended the ehective medium theory into Weyl semimetals with chiral 
anomaly. Our theory could give the large unsaturated magnetoresistance of Weyl semimet¬ 
als, like TaAs, NbAs, NbP etc. The power law at high magnetic held and turn-on magnetic 
held are studied and in qualitatively agreement with experimental data. Furthermore, we 
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investigated the angle-dependence MR. Theoretically, there is a transition from postive MR 
to negative MR when the chiral anomaly is not negligible to compare with the large pos¬ 
itive magnetoresistance background. The extended EMT could give a good explanation to 
magnetotransport of current experiments on Weyl semimetals candidates TaAs, NbAs, NbP 
etc. 
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FIG. 1. Structure and magnetoresistance of TaAs family, (a)Body-centred tetragonal structure of 
Ta(Nb)As(P). Ta(Nb) and As(P) atoms are represented by brown and green atoms, (b) The first 
Brillouin zone (BZ) in the momentum space reprinted from Ref.^. Red and blue dots stand for Weyl 
nodes with opposite chirality, (c) Side view of BZ for Fermi surface. The Weyl nodes in the kz=0 
plane form hole pockets, while the others form electron pockets symmetrically, (d)Experimental 
data on magnetoresistance of TaAs faimly at low temperature, TaAs (red,blue) from Ref.^, 
NbP(magenta) from Ref.^, NbAs(Olivia) from Ref.^, and WTe 2 (black) from Ref.^^. (d)Power 
law of the unsaturated magnetoresistance, the square dependence from carrier compensation(black 
dotted), the exactly linear dependence from Abrikosovs quantum limit (red solid), the quasi-linear 
dependence from duality theorem on 2D electron systems (Olivia dashed), and saturated MR of 
classical prediction (blue short-dashed). The shadow region are well predicted by the effective- 
medium theory. Inset shows the noramlized MR of TaAs family are quasi-linear above at IT, while 
WTe2 are perfectly square dependence. 
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FIG. 2. The magnetoresistance (MR) on perfect n-p compensation condition, (a) MR as a func¬ 
tion of magnetic field at different mobilities from 5 x to 10®(b) The 

corresponding Hall resistance pxV are all zeros at fixed n=p.(c)The electron conductivity and hole 
conductivity are equal, and the electron mobility is fixed at 10^cm^V~^s~^. MR as a function of 
hole mobility at different magnetic field lOT, SOT and GOT. 
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FIG. 3. The magnetotransport of Weyl semimetal at n greater than p or fin greater than //p.(a)-(c) 
at fixed fLn,p = crn?‘V~^s~^. (a)The longitudinal resistance pxx as a function of magnetic field 
at different carrier ratios p/n. (b) the corresponding transverse resistance pxy as a function of 
magnetic field. (c)The ”turn-on” magnetic field (black line) and the power law at high magnetic 
field (red line) as a function of the carrier ratio p/n from 0.1 to 1. (d)-(f) at fixed n = p and 
Pn = cm?V~^s~^. (d)The longitudinal resistance pxx as a function of magnetic field at different 
mobility ratios Pp/pn- (e) the corresponding transverse resistance pxy as a function of magnetic 
field. (f)The ”turn-on” magnetic field (black line) and the power law at high magnetic field (red 
line) as a function of the mobility ratio p/n from 0.1 to 1. 


11 











Ah 


AA'" AA 


CD 




(T) 



(T) 


FIG. 4. The angle-dependence magnetoresistance from positve MR to negative MR. (a)The chiral 
energy difference induced by chiral anomaly within parellel magnetic field H and electric field E. 
(b)The MR (points) as a function of magnetic field at different angle from 0 to tt/2 and fixed 
mobility fjLn,p = 10®The dashed lines show the positive MR without chiral anomaly as 
a reference, (c) The density plot of MR as a function of magnetic field and angle. The yellow line 
shows a transition between PMR to NMR at 5 = 0.3. (d) The contour plot of magnetoresistance 
transition at different parameter 5. The NMR region enlarges with 6 increasing. 
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